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(54) Post-processing method for reducing artifacts in block-coded digital images, and post- 
processing device for actuating such method 



(57) A post-processing method for reducing arti- 
facts in block-coded digital images, characterized by 
providing for: dividing an input image (In) into a plurality 
of image blocks (IB); for each image block (IB), estimat- 
ing global features (GF) of said image block providing 
information on an average content of image edges 
along the horizontal and vertical directions of said 
image block (IB); for each pixel (Px) of an image block 
(IB) under examination, estimating local features (LF) 



for said pixel (Px) providing information on the content 
of image edges along the horizontal and vertical direc- 
tions of an image area near said pixel (Px); modifying 
the value of said pixel (Px) according to both said global 
features (GF) of the image block (IB) to which said pixel 
(Px) belongs and said local features (LF) of the image 
area near said pixel (Px). 
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Description 

The present invention relates digital image coding. More precisely, the invention relates to a post-processing 
method for reducing artifacts in block-coded digital images, and to a post-processing device suitable for actuating such 
5 a method. 

With the diffusion of digital communication systems, digital images are more and more used. This has led to the 
diffusion of still and video cameras with digital acquisition and processing capability. 

In order to better exploit storage means and transmission bandwidth, digital image compression standards have 
been developed, such as JPEG for still images, and MPEG-1 and MPEG-2 for digital television image sequences. 
10 The above-referred compression standards provides for block-coding based on Discrete Cosine Transform (DCT). 
A digital image is divided into blocks of pixels, and each block is encoded independently from the others. DCT coeffi- 
cients for the pixels of each block are evaluated and a quantization matrix is applied to the DCT coefficients to reduce 
the information to be stored or transmitted. When the image is to be displayed, it must be decoded in advance. 

Due to the quantization process, these image compression methods are lossy, i.e. they cause a loss of information 
75 in the decoded image with respect to the original image. The decoded image can thus present noticeable degradation, 
mainly consisting of two kinds of artifacts known in the art under the names of "grid noise* and "staircase noise". 

In order to reduce the image degradation, post-processing methods of the decoded image have been proposed 
which allow for attenuating grid noise and staircase noise. 

In view of the state of the art described, it is an object of the present invention to provide a new post-processing 
20 method for reducing artifacts in block-coded digital images. 

According to the present invention, such object is attained by means of a post-processing method for reducing arti- 
facts in block-coded digital images, characterized by providing for: 

a) dividing an input image into a plurality of image blocks; 
25 b) for each image block, estimating global features of said image block providing information on an average content 

of image edges along the horizontal and vertical directions of said image block; 

c) for each pixel of an image block under examination, estimating local features for said pixel providing information 
on the content of image edges along the horizontal and vertical directions of an image area around said pixel; 

d) modifying the value of said pixel according to both said global features of the image block to which said pixel 
30 belongs and said local features of the image area around said pixel. 

Also according to the invention, there is also provided a post-processing device for actuating said method, charac- 
terized by comprising: 

35 • first means supplied with an input image for estimating global features of an image block under examination, said 
global features providing information on an average content of image edges along the horizontal and vertical direc- 
tions of said image block; 

second means supplied with said input image for estimating local features for each pixel of the image block under 
examination, said local features providing information on the content of image edges along the horizontal and ver- 
40 tical directions of an image area around said pixel; 

third means supplied with said global features and said local features for modifying the value of said pixel according 
to both said global features and said local features. 

The features and advantages of the present invention will be made apparent from the following detailed description 
45 of an embodiment thereof, illustrated as a non-limiting example in the annexed drawings, wherein: 

Figure 1 is a schematic block diagram illustrating the principle of operation of a method according to the present 
invention; 

Figure 2 shows a digital image divided into image blocks; 
so Figure 3 shows in detail an image block into which the digital image of Figure 2 is divided; 

Figure 4 shows an array of pixels of the image block of Figure 3; 

Figure 5 shows an image sub-block of the image block of Figure 3 used for evaluating global features of the image 
block; 

Figure 6 shows an horizontal processing window used for evaluating local features in the horizontal direction for a 
55 generic pixel of the image block; 

Figure 7 shows a vertical processing window used for evaluating local features in the vertical direction for said 
generic pixel; 

Figures 8 and 9 shows two membership functions used to perform a fuzzy computation; 
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Figure 10 is a block diagram of a device according to the present invention; 
Figure 1 1 shows the structure of two blocks of the device of Figure 10; and 
Figure 12 is a block diagram of other two blocks of the device of Figure 10. 

With reference to Figure 1 , there is shown a block diagram illustrating the principle of operation of the post-process- 
ing method according to the present invention. An input decoded compressed digital image In is supplied to a Feature 
Extraction block FE. Block FE provides for analyzing the image to evaluate global and local features thereof. The global 
and local features, respectively GF and LF, of the image In evaluated by block FE are supplied as inputs to a Fuzzy 
Process block FUZZY which, according to fuzzy rules, determines parameters FA suitable for determining the kind of 
filtering to be performed, in accordance to the global and local features GF and FL of the image. The parameters FA 
calculated by block FUZZY are supplied to a Filter Composition block FC which according to said parameters FA deter- 
mines the type of filtering to performed out of a set of predefined filters (block FS). Filter parameters FP determined by 
block FC are then supplied to a Processing block PROC, also supplied directly with the input image In, which performs 
the filtering of the input image In according to the filter parameters FP to provide a post-processed output image Out. 

It appears that the kind of filtering to be performed on the decoded input image In is chosen after an estimation of 
the global and local features of the decoded input image. For image areas near grid noise and near an edge, a low-pass 
filtering is performed, to reduce both staircase noise and grid noise. For areas containing fine details (image edges and 
texture), no filtering is performed. Thus, the method according to the present invention provides for performing a non- 
linear adaptive filtering on the pixels of the decoded image. 

The principle of operation outlined above will be now described in detail. 

As shown in Figure 2, the input image In is partitioned into image blocks IB, each containing an equal number of 
pixels. A typical dimension of the blocks is 8*8 pixels (Figure 3), but this is not however to be intended as a limitation, 
since other block dimensions are suitable. 

The image blocks IB of the input image In are scanned line by line starting from the top-left block to the bottom-right 
one. For each image block IB. the Feature Extraction block FE in Figure 1 determines the global and local features GF 
and LF. 

Global features of the image block IB under examination are determined by applying horizontal and vertical Sobel 
operators: 



horizontal Sobel operator (Hsob): 



rh11 h12 h13' 
h21 h22 h33 
h31 h32 h33. 



vertical Sobel operator (Vsob) : 



V11 v12 v13" 
v21 v22 v23 
v31 v32 v33. 



to each pixel belonging to an image sub-block internal to the image block IB. For example, the following Sobel opera- 
tors: 



Hsob: 



-1 0 V 
-2 0 2 
-1 0 1. 



Vsob: 



-1 -2 -V 

0 0 0 

1 2 1 



are applied to a 6*6 pixel image sub-block ISB (gray area in Figure 5). As shown in Figure 4, for each pixel Pi of the 
image sub-block ISB a 3*3 array of neighboring pixels M centered in pixel Pi is considered, and the values of the pixels 
of said array M are multiplied by the coefficients of the horizontal and vertical Sobel operators, to obtain: 

Hsob = (P3+P8+2*P5) - (P1+P6+2*P4), 
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Vsob = (P6+P8+2*P7) - (P1+P3+2*P2). 
where PI -P4 and P5-P8 are the values of the pixels (gray levels) 



taanlv^c^^ 

taJI* Ut V3lUeS ° f horizontal Sobel operators calculated for the pixels of image sub-block ISB are accumu- 
lated to obtain an accumulated value Acc(Hsob). and the output values of the vertical Sobel operators 
p.xels of image sub-block ISB are accumulated to obtain an accumulated value Acc(vsob). aS^T5^.^5T 
frequency content ,n the horizonta. direction (vertical edges) of the image block «B. Acc(Vsob) g JesThe K^5£, 
content ,n the vertcal d,rection (horizontal edges) of image block IB. Thus. Acc(Hsob) and ^i^nSSS^ 
vide the degree of "edgeness" of the image block under examination in the vertcal and horizonta? S , o e 
no ed that in order to evaluate the g.obal features GF of the image block IB under examination ZthS 
om.s block are considered (by applying 3*3 Sobel operators to the 6*6 image sub-block .SB, ittnSceiry °o take 
into consideration pixels belonging to neighboring image blocks) noi necessary to taKe 

a ,? 0 K a x' f !f tUfeS GF °' the image bl0Ck under exami "ation can be formed by the accumulated values Accf Hsob) and 

t0 ° bta,n ^ nUmbSr ° f ed96S in the hori - tal and erections in the 

hi ^ f ! atUr6S LF ° f * he ima9e b ' OCk IB are estimated according to the following method. All the pixels of the imaae 
block IB under examination are scanned .ine by line starting from the top-left pixel down to the bottom one T 0 ™?h 
pixe a honzonta processing window is applied: a prescribed number of pixels respectively preceding anS Sowing the 
pixel under cons^eration and belonging to the same image line of the pixel under consideration is ^^7^ ta- 
ble horizontal processing window HPW is shown in Figure 6, which is a horizonta. 1 *5 processing window for a Zen 
pixel, the two preceding p,xe.s Pa. Pb and the two following pixels Pc. Pd belonging to the same line a^nsid erST 
Figure 6 there is shown by way of example the horizontal processing window HPW associated to he rZ P ?x2 Px of the 
■mage block. It should be noted that not on.y the pixels of the image block IB under examination are cvSrteZ bufa so 
pjxds belonging to neighboring image blocks; this is for example the case of the first. sel . s ele i andligW ptel 
of each line of pixels of the image block IB under examination 

Similarly a vertical processing window is applied to each pixel of the image block IB The vertical orocessino win 

tZZ ^^aS^rf d '? ,0n PX -r 3 nUmb6r * ** be ' 0nging ^ SiSSTi 

ana preceding and following the pixel under consideration; for example, as shown in Figure 7 the vertical orocessinn 

window VPW can have dimensions identical to the horizonta. processing window HPW (5*1 )',^!SS£S^S- 
els Pe Pf preceding p.xel Px and two pixels Pg. Ph following pixel Px in the vertical direction P 
dow VPw!n nl n° be ' ?<* Pr&,k3US,y mentioned is ^ applied to each pixel of the vertical processing win- 

Slunl^ 

^SltiSS, high-frequency content in the vertical direction of an image region around the pixel 

Acc^^TZ? F ^ the ima9S b ' 0Ck ' B Und6r examination ( i e - »• two accumulated values Acc(Hsob) and 
nSl * ] T alternat,v e. »he average value of Acc(Hsob) and Acc(Vsob)) and the local features LF for the pixel 

^^ss^^is^ b f k / lhe te r values HS1 - HS5 and vsi - vs5) are then ««*- to"he he Fu px z 

process block FUZZY. Th.s block prov.des for evaluating the degree of membership of a generic value HSi and Vsi 
(.=1.-5) to two fuzzy sets "Small" and "Big"; these degree of membership can be evaluated Sy app^g to HS VSi the 

X^SEZZTZ in H R9UreS 8 a ' n *~* <i9UreS ' ™ ^ ™ are ^SeZl on *e ™Z 
1™! 21 T ? UndSr examination - j e- on the accumulated values Acc(Hsob) and Acc(Vsob) or on the 

case, i m ana i n2 are the same for Hsi and Vsi values 

sets -SS »R Vi - 9 35 anteced , en,S thS deQreeS ° f membershi P ° f output values HSi and VSi to the two fuzzy 
sets Small and Big are then evaluated. This means that 32 rules are to be evaluated for both the horizontal and vll 

S^Mnthis ::t a " ; hose ,u2z , y ru,es ha ™ 9 the same - * ™ £^ 

?h^.I ? T complexit y is reduced - and a total of nine rules for each direction have to be evaluate? 
The following fuzzy rules are applied to the five values HS1-HS5 associated to the horizontal direction 

1. If HSI is Small and HS2 is Small and HS3 is Small and HS4 is Small andHSS fs Small, then a1 is Big; 

2. IfHSI js Small andHSZ ls Small and HS3 is Small and HS4 is Small and HS5 is Big, then a2 is Big; 
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3. // HS1 is Small and HS2 is Small and HS3 is Small andHS4 is Big and HS5 is Small, /ten a3 is Big; 

4. //HS1 is Small and HS2 is Small and HS3 is Small andHSA is Big and HS5 /s Big, fnen <x4 /s Big; 

5. /f HS1 is Small and HS2 is Big and HS3 is Small and HS4 is Small and HS5 is Small, fnen a5 /s Big; 

6. //HS1 is Big and HS2 is Small and H$3 is Small and HS4 is Small and HS5 is Small, /nen a6 /s Big; 

7. // HS1 is Big and HS2 is Small and HS3 /s Small and HS4 /s Small and HS5 is Big, fnen a7 is Big; 

8. // HS1 is Big and HS2 /s Big and HS3 is Small and HS4 is Small and HS5 is Small, fnen a8 is Big. 

The activation level of each rule depends on the degree of memberships of the pattern of output values HSi of the 
horizontal Sobel operator applied to the five pixels of the horizontal processing window HPW; the degree of member- 
ships depend in turn on the global features GF of the image block to which the pixel under examination belongs. The 
activation level of the else (ninth) rule is computed as aelse = (1 - aave) , where aave is the average activation degree 
of fuzzy rules 1 to 8. al to a8 and aelse, and a similar set of nine activation degrees for the fuzzy rules applied to values 
VST -VS5) form the output FA of the fuzzy process block FUZZY in Figure 1 . 

Each one of the above-listed rules is associated to a respective set of predefined filter parameters, which are stored 
as a look-up table in block FS of Figure 1 . Suitable predefined filter parameter sets are for example: 

Rule 1: (c1 1=1.0. d 2=1.0, c1 3=1.0, c1 4=1.0, d 5=1.0) if the pixel under examination lies outside the image sub- 
block ISB, and (c11=0.0, c12=1.0, c13=1.0, c14=1.0, cl5=0.0) if the pixel under examination lies inside the image 
sub-block ISB; 

Rule 2: (c21=0.5, c22=1.0, c23=1.0, c24=1.0, c25=0.0); 
Rule 3: (c31=0.5, c32=1.0, c33=1.0, c34=0.0, c35=0.0); 
Rule 4: (c41=0.5. c42=1.0. c43=1.0, c44=0.0, c45=0.0); 
Rule 5: (c51=0.0, c52=0.0, c53=1.0, c54=1.0, c55=0.5); 
Rule 6: (c61=0.0, c62=1.0, c63=1.0. c64=1.0, c65=0.5); 
Rule 7: (c71 =0.0, c72=1 .-0. c73=1 .0, c74=1 .0, c75=0.0); 
Rule 8: (c81=0.0, c82=0.0. c83=1.0. c84=1.0, c85=0.5); 
Else rule: (c91=0.0, c92=0.0. c93=1 .0. c94=0.0. c95=0.0). 

The parameters FP of the filter to be applied to the five pixels of the horizontal processing window HPW are calcu- 
lated as a weighted average of the nine filters described above, with weight factors formed by the activation degrees al 
to a8 and aelse of the respective fuzzy rules. 

Assuming that ai is the activation degree of the i-th fuzzy rule (i=1 ..9), the.ninth fuzzy rule being the else fuzzy rule 
(a9 = aelse ), and cij are the coefficients of the i-th filter (i=1 ..9, j=1 ..5), the weight factor applied to the i-th filter, asso- 
ciated to the i-th fuzzy rule is: 

Fi = ai * Cij 

and the coefficients Hj of the final horizontal filter to be applied to the pixels of the horizontal processing window HPW 
are given by: 



where N is a normalization factor. 

The horizontally-filtered value Px of the pixel Px under examination (at the center of the horizontal processing win- 
dow) is then calculated as a weighted average of the values of the pixels Pa, Pb, Px, Pc and Pd belonging to the hori- 
zontal processing window HPW, with weight factors formed by the coefficients Hj: 

Px = H1 *Pa+ H2*Pb+ H3*Px+ H4*Pc+ H5*Pd. 

Similar calculations are performed for the vertical direction, starting from the output values VS1 -VS5 of the vertical 
Sobel operators applied to the pixels Pe. Pf. Px, Pg and Ph in the vertical processing window VPW. The coefficients Vj 
(j=1..5) of the filter for the vertical direction are calculated in a way completely similar to that used for determining the 
coefficients Hj: 
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i=1 



N 



10 
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vertical direction. The coefficSnte vj i areTen ^SlSSSC^J^T f^*" 1 t0 the ^ fU2zy rule for the 
a weighted average of the same The WerS vafue «Z nil* p V processin 9 wind °" VPW to calculate 

vertical direction, is proved at th^u^ 

o^TeZZZZ EE S inte" JS^F** ^ ^ filt " «^ V3 can be the value Px 
U ^ Ine p,xels ,n the n ° n zontal processing window HPW the horizontal filter HJ (j=i.. 5 ) : ~~ 

Out = Vl*Pe + V2*Pf + V3*Px + V4*Pg + V5*Ph. 
Alternatively, it is possible to eva.uate first the vertically-filtered value £* of the pixel under examination: 

Px = VI *Pe + V2*Pf + V3*Px + V4*Pg + V5*Ph. 

Out = H1 *Pa + H2*Pb + H3*Px + H4*Pc + H5*Pd. 

processed is dMded in, and a s^ond 2? ? ,? ° 01 * e ima9e """^ IB «» im =9» *> °» PM- 

two line memories LM1 andLM2 sr^ies "b£F f " S " - * ^ 8,50 5UPP ' ieS a <" 

6-6 image sub-Mock ISB inside ««S?matUorS IB t£ ? ^ ? "S S ° bel operators tor ,n « °> »» 
vsob of fh. hctan* and -^Ss^S^^u^JS^^ 4 the or^uts Hsob and 

r^a^^^ 

~ ^ sJ^^a^^^^^ 

in palaw S'«°e iiToT]!^^^,!^ ?' 5? W ""^ LM3-L M ,o. Block 2 is supp»ed 

pixei inside the vertica. P'oo^Z^^^l^^i^T V*" ** " M ^ ° perator to each 
ferred providing for calculating five venS Sab* «^ 1 ? J * 6 memories - a P^Mel approach is pre- 
VS1 -VS5 are supplied to a vert S C L^Sf m """Si ^ ° UtPUtS ° f the ,ive vertical Sobel 
7 and the output MOUT of the memory MEM of Nock 2S S^f . ^ ° ° f ,he pixel delav block 

currently processed by block 2 ™ the a^ZZ va Z Z^ * ^ ° F ° f *" im39e b '° Ck ' B 

Acc(Hsob). The vertical fuzzy f | ter block 9 e^tS £l ! Acc < Vs °b) or alternately, the average of Acc(Vsob) and 
"Small" and "Big" taking into accou^ ° f "alu- v S1-VS5 to the fuzzy sets 

fuzzy rules for the vertical direction. iS^SliSS^^^ S^L^Tu T* ^ ° f the nine 

S-stheveWe^reS^^^ 
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plies directly a horizontal fuzzy filter block 10; the output Ed of block 9 also supplies a cascade of four pixel delays D 
whose outputs respectively form the vertically-filtered values Eg, E& Eh, Ea of the pixels Pc, Px. Pb, Pa in the horizontal 
processing window HPW and supply the horizontal fuzzy filter block 10. 

In parallel to blocks 8 and 9, the outputs L7-L9 of the pixel delay block 7 supply a block 1 1 which applies the hori- 
zontal Sobel operators to the pixels inside the horizontal processing window HPW. Differently from the vertical sobel 
operators, only one horizontal sobel operator is calculated at a time; a compensation delay block 12 introduces a delay 
for compensating the processing delay of the vertical fuzzy filter block 9. The output of block 12 t forming the output of 
the horizontal Sobel operator HS5 applied to pixel Pd of the horizontal processing window in Figure 6, supplies the hor- 
izontal fuzzy filter block 10 and a cascade of four pixel delays D, the outputs thereof forming the values HS4. HS3, HS2 
and HS1 and suppling the horizontal fuzzy filter block 10. The horizontal fuzzy filter block 10, which is also supplied by 
the output MOUT of the memory block MEM in block 1 providing the value Acc(Hsob) (or alternatively the average of 
values Acc(Hsob) and Acc(Vsob)), evaluates the degree of membership of values HS1-HS5 to the fuzzy sets "Small" 
and "Big" according to the value of the global features GF provided by MOUT, evaluates the activation levels of the nine 
fuzzy rules described above for the filtering in the horizontal direction, calculates the coefficients Hj of the horizontal fil- 
ter and applies the parameters Hj to the vertically-filtered values E§, Efe. Ex.. E£, Ed of the pixels Pa. Pb, Px, Pc. Pd in 
the horizontal processing window HPW to obtain the horizontally- and vertically-filtered value Out of the pixel Px under 
examination, 

A control circuit CTRL controls the operation of blocks 1 , 2 and the line memories LM1-LM10. 

Figure 1 1 shows the structure of the vertical and horizontal Sobel operator blocks 8 and 1 1 of Figure 10. They are 
composed in a straightforward way by adders. 

Figure 12 shows the structure of both the vertical fuzzy filter block S and the horizontal fuzzy filter block 10. X1 -X5 
are the vertical or, respectively, horizontal Sobel operator outputs VS1-VS5 and HS1-HS5. X1-X5 are supplied to a 
fuzzy rule evaluation block 1 3 which evaluates the activation degrees pi -p9 of the nine fuzzy rules for the vertical direc- 
tion or, respectively, the activation degrees a1-a9 of the nine fuzzy rules for the horizontal direction. The activation 
degrees evaluated by block 13 are supplied to a look-up table of respective predefined filter parameters F1-F9 (forming 
block FS in Figure 1). and the outputs of the look-up table, i.e. the predefined filter parameters cij multiplied by the acti- 
vation degree of the respective fuzzy rule, are supplied to a filter composition block 14 which calculates the coefficients 
V1-V5 or, respectively. H1-H5. of the vertical or, respectively, horizontal filter. Said coefficients are then supplied to a 
processing block 15 which is also supplied with the pixel values PXS (L6-L10 or, respectively, P§, P^, Px^, Pc_, Ed in Fig- 
ure 10); block 1 5 applies the filter coefficients to the pixel values to obtain the filtered value of the pixel under examina- 
tion Px. 

Claims 

1 . Post-processing method for reducing artifacts in block-coded digital images, characterized by providing for: 

a) dividing an input image (In) into a plurality of image blocks (IB); 

b) for each image block (IB), estimating global features (GF) of said image block providing information on an 
average content of image edges along the horizontal and vertical directions of said image block (IB); 

c) for each pixel (Px) of an image block (IB) under examination, estimating local features (LF) for said pixel (Px) 
providing information on the content of image edges along the horizontal and vertical directions of an image 
area near said pixel (Px); 

d) modifying the value of said pixel (Px) according to both said global features (GF) of the image block (IB) to 
which said pixel (Px) belongs and said local features (LF) of the image area near said pixel (Px). 

2. Method according to claim 1 , characterized in that step d) provides for: 

d1) defining a set of predefined local features; 

d2) determining degrees of coincidence (al -a9,pi-p9)of said local features (LF) of the image area around said 
pixel (Px) with each predefined local features of said set, said degrees of coincidence (a1-a9,pi-p9) depend- 
ing on said global features (GF) of the image block (IB) to which said pixel (Px) belongs; 
d3) making the value of said pixel (Px) equal to a weighted average (Ex.ExJ of the value of said pixel (Px) and 
of the values of neighboring pixels (Pa-Pd,Pe-Pf). with weight factors (H1-H5.V1-V5) depending on said 
degrees of coincidence of said local features (LF) with each of said predefined local features. 

3. Method according to claim 2. characterized in that said determining the degrees of coincidence (a1-a9,pi-p9) in 
step d2) provides for performing a fuzzy calculation. 
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said image sub-block (ISB) to obtain an accumui^H IS f ?L 6 ' ( Hs °») applied to each pixel of 

outputs of the vertical SoM operSors^b) ap^ 2 J TrZi TT*' **" ° Pera, ° rS and a3 ng 

mulated output of vertica. Sobe. operator^?^ P *' ° f ^ ,ma96 SUb - b,0ck < ISB > * aw? 

r r ir^r^r r ? the — b, ~ k « — 

(Acc(Hsob).A CC (Vsob)). ^ 3UtS " the h ° r| zontal and vertical Sobel operators 

^^ZVT^^^JZl'T, T - — «°> — -* 

(Acc(Hsob).Acc(Vsob)). accumulated outputs of the horizontal and vertical Sobel operators 

thepixel (Px); 9 ,nQ to a ^e ,mage l.neasthep,xel(Px)and preceding and following 

c2) applying said horizontal Sobel ooeratnr tuenM «~ „ u • . ,~ 

SSSS^SS ^e^^^ -ca, processing window 

set "Big" evaluating activation degrees S^oi^^T , ^ ^ SBt " Sma "" and to a ** 
one predefined horizontal pattern of horizontal ^Toa^ZTi: ^ ^ ^ 006 associ ^d with at least 
each vertica. Sobe. operator output (VS1-VS5 of SlSSiS^ 'T,™ 9 dWS ° f •""*—*> of 
fuzzy set "Big" and evaluating activation <^ZwE^Z?F 7 ** to 3 Second 

a. least one predefined vertical pattern of vertical L^J a ^Xutl * "* aSS ° Ciated with 

^"pS^^ of membership of the horizont, 

determining a first and a second trLbensi^Z^l^ „ * S6tS " Sma "" and " Bi 9" P™'^ for 

ing degrees of membership of the JSSi^^^^T^^ ^ ^ Said dete ™" 
toy 1*. "Smar and "Big" provides for 0***^*!^^^ ? V6rtical P 3 "^ to «W second 
global features (GF). '^mining a third and fourth membership functions depending on said 
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set being associated with a respective one of said groups of predefined horizontal weight factors, and each fuzzy 
rule of said second set being associated with a respective one of said groups of predefined vertical weight factors. 

14. Method according to claim 13. characterized in that said weight factors (H1-H5.V1-V5) comprise horizontal weight 
5 factors (H1-H5) and vertical weight factors (V1-V5), said horizontal weight factors (H1-H5) being determined by 

making a weighted average of the predefined horizontal weight factors with weight coefficients being formed by the 
activation degrees (a1-a9) of the fuzzy rules of the first set, and said vertical weight factors (V1-V5) being deter- 
mined by making a weighted average of the predefined vertical weight factors with weight coefficients formed by 
the activation degrees (p1 -p9) of the fuzzy rules of the second set. 

w 

15. Method according to claim 14, characterized in that the value of the pixel (Px) under examination is modified by 
applying the horizontal weight factors (H1-H5) to the values of the pixels (Pa-Pd,Px) in the horizontal processing 
window (HP W) and applying the vertical weight factors (V1 -V5) to the values of the pixels (Pe-Ph,Px) in the vertical 
processing window (VPW). 

15 

16. Post-processing device for reducing artifacts in block-coded digital images, characterized by comprising: 

first means (1) supplied with an input image (In) for estimating global features (MOUT) of an image block (IB) 
under examination, said global features providing information on an average content of image edges along the 

20 horizontal and vertical directions of said image block (IB); 

second means (7,8,1 1.12, D) supplied with said input image (In) for estimating local features (VS1-VS5.HS1- 
HS5) for each pixel (Px) of the image block (IB) under examination, said local features providing information on 
the content of image edges along the horizontal and vertical directions of an image area around said pixel (Px); 
third means (9,D,10) supplied with said global features (MOUT) and said local features (VS1-VS5.HS1-HS5) 

25 for modifying the value of said pixel (Px) according to both said global features (MOUT) and said local features 

(VS1-VS5.HS1-HS5). 

17. Device according to claim 16, characterized in that said first means (1) comprises means (3.4) for evaluating hori- 
zontal and vertical Sobel operator outputs (Hsob.Vsob) of horizontal and vertical Sobel operators applied to pixels 

30 of an image sub-block (ISB) internal to said image block (IB) under examination, and accumulator means (5) for 
accumulating the horizontal Sobel operator outputs (Hsob) and the vertical Sobel operator outputs (Vsob) for each 
pixel of the image sub-block (ISB). 

18. Device according to claim 17, characterized in that said second means (7,8,1 1,12, D) comprises fourth means (7,8) 
35 for evaluating vertical Sobel operator outputs (VS1-VS5) of vertical Sobel operators (Vsob) applied to said pixel 

(Px) and to vertically-neighboring pixels (Pe.Pf.Pg.Ph) preceding and following the pixel (Px) in the vertical direc- 
tion, and fifth means (11,12,D) for evaluating horizontal Sobel operator outputs (HS1-HS5) of horizontal Sobel 
operators (Hsob) applied to said pixel (Px) and to horizontally-neighboring pixels (Pa.Pb.Pc.Pd) preceding and fol- 
lowing the pixel (Px) in the horizontal direction. 

40 

19. Device according to claim 18. characterized in that said third means (9.D.10) comprises vertical filter means (9.D) 
supplied with said vertical Sobel operator outputs (VS1 -VS5) and said global features (MOUT) for calculating a ver- 
tically-filtered value ( Px) of said pixel (Px) depending on said global features (MOUT), said vertical Sobel operator 
outputs (VS1-VS5) and the values of said vertically-neighboring pixels (Pe.Pf.Pg.Ph), and horizontal filter means 

45 . (10) supplied with said global features (MOUT), said horizontal Sobel operator outputs (HS1-HS5) and said verti- 
cally-filtered value ( Px) of said pixel (Px) for calculating a horizontally-filtered value of said pixel (Px) depending on 
said global features (MOUT), said horizontal Sobel operator outputs (HS1-HS5) and the vertically-filtered values 
(Ei'SkE&Efi.Ed) of said P ixel ( Px ) and horizontally-neighboring pixels (Pa.Pb.Pc.Pd). 

so 20. Device according to claim 1 9, characterized in that said vertical filter means (9, D) comprises first fuzzy computation 
means (13) supplied with said vertical Sobel operator outputs (VS1-VS5) and said global features (MOUT) for eval- 
uating degrees of coincidence (pi -p9) of the pattern of vertical Sobel operator outputs (VS1-VS5) with a set of pre- 
defined patterns of vertical Sobel operator outputs, said degrees of coincidence (P1-p9) depending on said global 
features (MOUT). a look-up table of predefined vertical filter coefficients (F1-F9) each one associated to at least 

55 one respective predefined pattern of vertical Sobel operator outputs, a vertical filter coefficients composition means 
(14) supplied with said predefined vertical filter coefficients for generating a set of vertical filter coefficients (V1 -V5) 
which are a weighted average of the predefined vertical filter coefficients with weight coefficients formed by said 
degrees of coincidence (p1-(J9), and a vertical filter (15) supplied with said vertical filter coefficients (V1-V5) and 



9 
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21. Device according to claim 20. characterized in that said horizontal filter means (1 0) comprises second f u2 zv r™, 

one assocated to at least one respecliv. predefined pattern ol horizomal Sobel opera,™ oSTs a taSSS 
coeftaents composhon means (14) supplied with said predefined horizontal filter coe«enS^ F«^™ ? 
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